Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN I^TeX style file v2.2) 



The fingerprint of binary intermediate mass black holes in 
globular clusters: supra- thermal stars and angular 
momentum alignment 



M. Mapelli\ M. Colpi^ A. Possenti^ & S. Sigurdsson^ 

^S.I.S.S.A., Via Beirut 2 - 4, 1-34014 Trieste, Italy; mapelliSsissa. it 

^ Dipartimento di Fisica G. Occhialini, Universita di Milano Bicocca, Piazza della Scienza 3. 1-20126 Milano, Italy; 
monica. colpiamib . inf n. it 

^INAF, Osservatorio Astronomico di Cagliari, Poggio dei Pini, Strada 54, 09012 Capoterra, Italy; possentiSca. astro . it 
'''Department of Astronomy and Astrophysics, The Pennsylvania State University, 525 Davey Lab, University Park, PA 16802; 
steinnSastro .psu. edu 



2 February 2008 



ABSTRACT 

We explore the signatures that a binary intermediate mass black hole (IMBH) imprints 
on the velocity and on the angular momentum of globular cluster stars. Simulating 
3-body encounters between a star and a binary IMBH, wc find that the binary IMBH 
generates a family of few hundreds of stars (~ 100-300) which remain bound to the 
globular cluster (GC) and have velocity significantly higher than the dispersion veloc- 
ity. For this reason we term them "supra-thermal" stars. We also notice that, after 
the interaction, a considerable fraction (55-70%) of stars tend to align their orbital 
angular momentum with the angular momentum of the binary IMBH, introducing an 
anisotropy in the angular momentum distribution of cluster stars. We simulate the 
dynamical evolution of these supra-thermal stars before thermalization, and find that 
these stars tend to cluster at a distance of few core radii from the GC center. We 
conclude that the detectability of such signatures appears problematic with present 
telescopes. 
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1 INTRODUCTION 

A number of different observations suggest tliat large black 
holes (BHs) may exist in nature, with masses between 
2OM0 — 10* M0. Heavier than the stellar-mass BHs born 
in core-collapse supernovae {ZMq — 2QMq\ Orosz 2002), 
these intermediate mass black holes (IMBHs) are expected 
to form from the direct collapse of very massive stars, or in 
dense stellar systems through complex dynamical processes. 
Filling the gap between stellar-mass black holes (BHs) and 
super-massive black holes (SMBHs), they are of crucial im- 
portance in estabUshing the potential physical link between 
these two classes. 

IMBHs may plausibly have formed in the early universe 
as remnants of the first generation of metal free stars (Abel, 
Bryan & Norman 2002; Heger et al. 2003; Schneider et al. 
2002). If this is true, IMBHs can participate the cosmic as- 
sembly of galaxies and be incorporated in larger and larger 
units, becoming seeds for the formation of SMBHs (Madau 
& Rees 2001; Volonteri, Haardt & Madau 2003). IMBH may 
still be forming in young dense star clusters vulnerable to 
unstable mass segregation, via collisions of massive stars 



(Portegies Zwart 2004; Portegies Zwart et al. 2004; Porte- 
gies Zwart & McMillan 2002; Giirkan et al. 2004; Freitag et 
al. 2005a, 2005b; see van der Marel 2004 for a review). Un- 
ambiguous detections of individual IMBHs do not exist yet, 
but there are observational hints in favor of their existence, 
from studies of ultra-luminous X-ray sources in nearby star- 
forming galaxies (Fabbiano 2004; Mushotzky 2004; Miller & 
Hamilton 2002). 

It has also been long suspected that globular clusters 
(GCs) may hide IMBHs in their cores (Shapiro 1977; Shapiro 
& Marchant 1978; Marchant & Shapiro 1979, 1980; Dun- 
can & Shapiro 1982; see Shapiro 1985 for a review), despite 
the fact that their formation root is poorly known: runaway 
mergers among the most massive stars, at the time of cluster 
formation, can lead to an IMBH, similarly to what has been 
conjectured to occur in young dense star clusters (Porte- 
gies Zwart et al. 2004; Fregeau et al. 2004). An alterna- 
tive pathway is based on the idea that IMBHs form later in 
the evolution of GCs through mergers of stellar-mass BHs. 
Segregating by dynamical friction in the core, these BHs 
are captured in binaries that form through dynamical en- 
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counters with stars and BHs. The picture is that hardening 
by subsequent interactions with BHs lead them to merge, 
emitting gravitational waves (Miller & Hamilton 2002). BHs 
more massive than stellar-meiss BHs are thus created. How- 
ever, the interactions that produce hardening also provide 
recoil, causing the ejection of BHs (Kulkarni, Hut & McMil- 
lan 1993; Sigurdsson & Hernquist 1993; Portegies Zwart & 
McMillan 2000), but the heaviest black holes may remain, if 
their mass exceeds a still uncertain value between ~ SOM© 
and ~ 3OOM0 (Miller & Hamihon 2002; Colpi, Mapelli & 
Possenti 2003; Gultekin, Miller & Hamilton 2004). 

At present, optical observations of GCs hint in favor 
of the existence of IMBHs. In particular Gcbhardt, Rich & 
Ho (2002) suggest the presence of a 2tliX x lO" M© IMBH, 
to explain the kinematics and the surface brightness profile 
of the globulax cluster Gl in M31. Gerssen et al. follow the 
same method to indicate the possible presence of a l.Tti'y ^ 
10'^ Mq IMBH in the galactic globular cluster M15 (Gerssen 
et al. 2002, 2003). These indications arc still controversial, 
as Baumgardt et al (2003a, 2003b) showed that both the 
measurements of M15 and those of Gl can be explained 
with a central concentration of compact objects (neutron 
stars and white dwarfs) instead of the presence of a massive 
black hole. However, a more recent paper by Gebhardt, Rich 
& Ho (2005) further supports the hypothesis of an IMBH in 
Gl. In addition, there has been some claim of the existence 
of rotation in the core of a few GCs (Gebhardt et al. 1995; 
Gcbhardt et a. 1997; Gebhardt et al. 2000) suggesting the 
presence of a steady source of angular momentum in their 
cores. 

Radio observations can also help in discovering con- 
centrations of under-luminous matter in the core of GCs, 
thanks to the presence of millisecond pulsars that can probe 
the underlying gravitational field of the cluster. Currently, 
there are more than 100 known millisecond pulsars in GCs 
(Possenti 2003; Ransom et al. 2005; Camilo & Rasio 2005) 
and some show large and negative period derivatives. This 
is an important peculiarity, because millisecond pulsars spin 
down intrinsically due to magnetic braking. Negative period 
derivatives come from the variable Doppler shift caused by 
the acceleration of the pulsar in the gravitational field of 
the cluster itself (Phinney 1993). When these negative pe- 
riod derivatives can be ascribed solely to the gravitational 
field of the cluster, these pulsars place a lower limit on the 
mass enclosed inside their projected distance from the glob- 
ular cluster center. 

Two pulsars in M15, at 1" from the cluster center were 
found (Phinney 1993) with negative values of their period 
derivatives, consistent with the mass distribution implied by 
the stellar kinematics. The recent discovery by D'Amico et 
al. (2002) of two accelerated pulsars in NGC 6752 at 6" and 
7" from the center has highlighted the presence of about 2 x 
10^ Mq of under- luminous matter in the core of the cluster 
(Ferraro et al. 2003a), making NGC 6752 a special target for 
the search of an IMBH. NGC 6752 is interesting also because 
it hosts two pulsars in its halo. PSR-A is a millisecond pulsar 
in a binary system with a white dwarf (Bassa et al. 2003; 
Ferraro et al. 2003b), located at 3.3 half mass radii away: it is 
the farthest pulsar ever observed in a cluster. Colpi, Possenti 
& Gualandris (2002), and Colpi, Mapelh & Possenti (2003, 
2004) have suggested that this pulsar has been propelled in 
the halo due to a dynamical interaction with a binary IMBH. 



In modeling the gravitational encounter, Colpi, Mapelli & 
Possenti (2003) found that a (50-200, 10) M© binary IMBH 
is the preferred target for imprinting the large kick to the 
pulsar. Colpi et al. (2005) shortly reviewed the dynamical 
effects that a binary IMBH would imprint on cluster stars. 

Considering all the recent hints provided by optical and 
radio observations, we explore, in this paper, an new way to 
unveil a binary IMBH in a GC, previously overlooked, that 
exploits the dynamical fingerprint left by a binary IMBH on 
cluster stars. In particular wo like to address a number of 
questions: (i) What signature does a binary IMBH imprint 
on cluster stars ? Are the stars heated during the scattering 
process still remaining bound to the cluster ? (ii) Is there 
direct transfer of angular momentum from the binary IMBH 
to the stars to produce some degree of alignment, given the 
large inertia of the BHs ? (iii) Are prograde or retrograde 
orbits equally scattered ? 

In this paper, we simulate 3-body encounters of cluster 
stars with a binary IMBH, studying the energy and angular 
momentum exchange, and reconstructing the trajectories of 
those stars that are scattered away from equilibrium by the 
binary IMBH. In Section 2 we outline the method used to 
simulate 3-body encounters between a binary IMBH and the 
cluster stars. In Section 3 we present our main results: the 
formation, in a GC hosting a binary IMBH, of a population 
of high velocity (bound) stars (that we call "suprarthermal 
stars"), which tend to align their orbital angular momentum 
with that of the binary IMBH. In Section 4 we investigate on 
the expected number of these supra-thermal stars. In par- 
ticular, we use an upgraded version of the code presented in 
Sigurdsson & Phinney (1995) to follow the dynamical evo- 
lution of these supra-thermal stars inside a GC model which 
reproduces the characteristics of NGC 6752. This study al- 
lows us to put some constraint on the detectability of this 
family of high velocity stars. Section 5 contains our conclu- 
sions. 



2 THE SIMULATIONS 

We simulated 3-body encounters involving a binary IMBH 
(composed of two BHs of mass Mi and M2) and a cluster 
star of mass m = 0.5 Mq. The dynamics of the encounter 
is followed solving the equations of motion with a rmmcri- 
cal code based on a Runge-Kutta fourth order integration 
scheme with adaptive stepsize and quality control (explained 
in Colpi, Mapelh & Possenti 2003, hereafter CMP). For the 
target binary IMBH we considered an interval of masses be- 
tween 60 M0 up to 210 Mq. These are the favored masses 
of the hypothetical binary IMBH in NGC 6752. The bi- 
nary has semi-major axis a of 1, 10, 100, and 1000 AU and 
eccentricity^ e = 0.7, which is nearly the average binary ec- 
centricity in statistical equilibrium (Hills 1975). The set of 
models is described in Table 1. 

The initial conditions, that are Monte Carlo generated 
as in CMP, are sampled using the prescriptions indicated in 
Hut & Bahcall (1983). In particular, the relative velocity w'" 

^ We have repeated select simulations for different values of the 
eccentricity and find only minor quantitative differences ( ^ 10%) 
in our results. 



The fingerprint of binary intermediate mass black holes 3 



is distributed homogeneously between 8.5 — 11 km s^^ (ac- 
cording to the value of stellar velocity dispersion measured 
in NGC 6752 which is our reference cluster; Dubath, Mey- 
lan & Mayor 1997), and the impact parameter b is drawn 
at random from a probability distribution uniform in 
and in a range going from to a truncation value bmax- 
The truncation value bmax is chosen by requiring that the 
simulations include all the encounters with non-negligible 
energetic exchange, i.e. with outgoing velocities of the in- 
teracting star significantly higher than the initial velocity 
(see Appendix A for a discussion of our choices of the im- 
pact parameter ranges). The three orientation angles and 
the phase of the binary arc generated as indicated in Ta- 
ble 1 and 2 of Hut & Bahcall (1983). We initiate (termi- 
nate) integration when the distance between the incoming 
(outcoming) star and the center of mass of the binary is 
comparable to the radius of gravitational influence of the 
IMBHs (ra ~ 2G(Mi + M2)/a'^, where a is the 1-D stel- 
lar dispersion velocity). At the start of each simulation, we 
place the binary BH and the single star on their respective 
hyperbolic trajectories. 

At the end of every scattering experiment we store the 
final binding energy E^'^ and angular momentum of 
the binary IMBH, together with the velocity at infinity (or 
post-encounter asymptotic velocity) u^'", defined as the ve- 
locity of the interacting star after the encounter, estrapo- 
lated at infinity, and the angular momentum of the 
interacting star. In this paper we want, in particular, to 
quantify the mean change in the absolute value J* of the 
vector J* and the extent of alignment of J* in the direction 
of Jbh- 



3 RESULTS 

3.1 Supra-thermal stars 

We find that a sizable number of high velocity stars is cre- 
ated, under certain conditions, that may highlight the pres- 
ence of an IMBH in the cluster. These stars gain kinetic 
energy and an excess velocity relative to the mean, remain- 
ing bound to the cluster. In Table 3 we give the fraction of 
bound stars, defined as those having a final velocity lower 
than the escape velocity from the core ~ 35-40 km s~^, 
and the fraction of high velocity or supra-thermal stars, de- 
fined as having a post-encounter asymptotic (or "at infin- 
ity") speed between 20 and 40 kms^^. Depending on the 
characteristics of the binary IMBH, these stars can account 
for ^ 50% of all stars that have experienced an encounter. 

Figure 1 shows the distribution of the velocity at infinity 
of stars scattering off a binary IMBH, for cases Dl, D2, D3 
and D4 (in Table 1). The shaded area indicates the strip 
of suprar-thermal stars. The mean values of the speed at 
infinity ii-'^*" and of the fractional binding energy exchange 
(ASbh/B^Th) = {{E^H - ESk)/ESk) are given in Table 2 
together with the dimensionless factor ^e, defined by: 



(AEm/ESii) = £.E [m/(Mi + M2)] . 



(1) 



As illustrated in Figure 1, dynamical encounters widen the 
velocity distribution of the stars. This effect is particularly 
severe when the binary is hard, having a higher binding en- 
ergy which becomes available in the interaction, (^e clusters 



around values between 0.6 and 4, and shows a maximum 
declining as the binary becomes extremely hard. This ex- 
plains why case Dl has a peaJj velocity at infinity smaller 
than case D2 corresponding to a less hard binary (see Ta- 
ble 2). We tried to fit the behavior of ^_b, as a function 
of the orbital separation a, with a parabola. We obtained 
= 6Log^(a) + cLog(a) -h d for 6=-1.2495, c=3.8615, 
c?=0.958; this fit is only an approximation, since we have 
few numerical points (Fig. 2). It is interesting to notice that 
for a = 1 AU is lower than for wider binaries. This docs 
not represent a violation of the first statement of Heggie's 
law (i.e. that hard binaries tend to become more energetic, 
Hcggie 1975; Hills 1990); in fact all our hard binaries be- 
come harder and harder. But a decreasing for a < 10 AU 
seems in partial disagreement with the second statement of 
Heggie's law, for which the hardening rate of the binaries is 
" approximately independent of their binding energy" (Heg- 
gie 1975). This difference is the result of a different approach 
with respect to that followed by Heggie. In fact, we are treat- 
ing very massive binaries with respect to the incoming star, 
and such binaries require very large maximum impact pa- 
rameters, nearly independently of their semi-major axis, to 
include in our simulations all the interactions which have 
significant energy exchange (i.e. ABBH/iSBH > 10~^), corre- 
sponding to post-encounter asymptotic velocities of the star 
in the supra-thermal range (see Appendix A for a complete 
discussion) . 

In Figure 1, case Dl and D3 give the highest rmmber of 
supra-thermal stars among the different runs. These supra- 
thermal stars axe the ones that have also acquired a siz- 
able fraction of the orbital angular momentum of the binary 
IMBH. Thus in scattering off the binary IMBH they are 
preferentially launched in a halo orbit, inside the cluster. 
Their detection (discussed in § 4) would be the sign of an 
IMBH hidden in the cluster. 



3.2 Angular Momentum Transfer and Alignment 

An interesting question to address is whether and how the 
orbital angular momentum of the binary IMBH couples to 
the star after a close gravitational encounter. The BHs can 
be sufficiently massive and the orbit be sufficiently wide that 
the orbital angular momentum of the binary exceeds that of 
the incoming star. In this case a direct transfer of orbital 
angular momentum can occur. The star coming close to the 
binary IMBH can be dragged into corotation, i.e. the star 
can emerge after the encounter with an angular momentum 
nearly aligned with the binary IMBH. 

If we denote with ^ = MiAh/iMi + M2) the reduced 
mass of the binary hosting the two BHs, the total angular 
momentum of the system is 



r m(Mi-HM2) 1 ^1 



(2) 



= Jbh + J* 

= /xV«G(Mi+M2)z + bu*- [-M,+M.+^J 

where z and z' are the unit vectors indicating respectively 
the directions of S^m and J*". Angular momentum transfer 
from the binary to the interacting star and partial alignment 
become important when ■Pgij ^ J™, i.e., when 



-^J{M^+M2)aG■>\. 



Our experiment confirms this fact. Figure 3 {solid lines) 
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Table 1. Initial Parameters. 







A/J' 


hq 


Semi- major axis (a) 


Range of impact parameters 


Number of simulations 


CASE 


Al 


50 


10 


1 


0-60 


5000 


CASE 


A2 


50 


10 


10 


0-60 


10000 


CASE 


A3 


50 


10 


100 


0-100 


10000 


CASE 


A4 


50 


10 


1000 


0-1000 


10000 


CASE 


Bl 


50 


50 


1 


0-100 


5000 


CASE 


B2 


50 


50 


10 


0-100 


5000 


CASE 


B3 


50 


50 


100 


0-100 


5000 


CASE 


B4 


50 


50 


1000 


0-1000 


5000 


CAbE 


CI 


100 


10 


1 


0-100 


5000 


CASE 


C2 


100 


10 


10 


0-100 


5000 


CASE 


C3 


100 


10 


100 


0-100 


5000 


CASE 


C4 


100 


10 


1000 


0-1000 


5000 


CASE 


Dl 


100 


50 


1 


0-100 


5000 


CASE 


D2 


100 


50 


10 


0-100 


5000 


CASE 


D3 


100 


50 


100 


0-100 


5000 


CASE 


D4 


100 


50 


1000 


0-1000 


5000 


CASE 


El 


200 


10 


1 


0-100 


5000 


CASE 


E2 


200 


10 


10 


0-100 


5000 


CASE 


E3 


200 


10 


100 


0-100 


5000 


CASE 


E4 


200 


10 


1000 


0-1000 


5000 



" In units of the solar mass. 
^ In Astronomical Units. 
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Figure 1. Post-encounter asymptotic velocity distributions of the 
cluster star in the case Dl {dotted- dashed line), D2 {dotted line), 
D3 {dashed line), D4 {solid line). The shaded area refers to the 
supra-thermal stars. On the y-axis the number of cases for each 
bin is normalized to the total number of resolved runs (Ntot). 



shows the post-encounter distribution of J* for cases D2, 
D3 and D4. It is compared with the distribution of J, of 
the incoming stars {dashed lines) to show that the binary 
IMBH transfers angular momentum to the stars widening 
the distribution of J, . Table 2 contains the averaged value 



Figure 2. = [{Mi + M2)/m\{AEBH/E^^} as a function of 
the orbital separation a for all the considered cases. Al, A2, A3 
and A4 are represented by open triangles; Bl, B2, B3 and B4 by 
crosses; CI, C2, C3 and C4 by open squares; Dl, D2, D3 and D4 
hj filled circles; El, E2, E3 and E4 by open circles. The solid line 
indicates a parabolic fit feLog^(a) 4- cLog(a) + d = ^e{o.), where 
6=-1.2495, c=3.8615, (i=0.958. 



of the fractional angular momentum increase (in modulus) 
(AJ,/J,'") for the complete series of runs. (AJ,/J™) is 
large and exceeds unity in correspondence to the runs where 
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Table 2. Final velocities and angular momentum exchanges.* 



.Jin 



AJ d 



CASE 


Al 


r:+95.0 

±^•0 Q 


0, 


.056 Q 


0, 


™r-+0.100 
■"^■^-0.100 


0, 


,oo86t°-.;;^?° 


1, 


.03 


CASE 


A2 


07 c+35.0 
— oU.U 


0, 


,630+?/24 


0, 


-17^+1.150 
■-^"^-0.850 


0, 


Q999+1.1032 


3 


.59 


CASE 


A3 


12.5+1°-° 

— lU.U 


1, 




0, 


„r,c. + 2.2()() 

■ '^^^-1.400 


0, 


rjr)77+1.0363 
■ \j ^ { i nmni 

— u.uiy i 


3, 


.32 


CASE 


A4 


7 c:+2.0 


0, 


-|c;c- + 3().33() 
' —0.228 


0, 


Q9t: + 0.35() 


0, 


^-,19+0.0309 
— 0.0199 


1, 


.34 


CASE 


Bl 


i2.5+^;^-» 


0, 


g.„ + 7.98() 
— U.Uo / 


0, 




0, 


QQO1+0.1778 

— \}.\}\}0\J 


0, 


.62 


CASE 


B2 


67.5+6°:° 


0, 


0-1 7+31.840 
. <J-L ' Q 34g 


0, 


„9t-+0.65 
■"^°-0.55 


0, 


■0168l°-.g?gi 


3, 


.36 


CASE 


B3 


17.5+25-° 
— 10. u 


1, 


994+145.935 

— i.40D 


0, 


O9r;+2.850 


0, 


— U.Ulo4 


3, 


.95 


CASE 


B4 




0, 


, 97+17.279 
■ f Q 259 


-0.02551°;^^° 


0, 


.0065 Q QQ73 


1, 


.30 


CASE 


CI 




0, 


n2n+5-''2® 
■"^"-0.020 


0, 


nntS+O.OlO 
■^"'^-0.005 


0, 


(,(,00+0.7152 

.UUOO_g QQ32 





.73 


CASE 


C2 


or, c: + 40.0 


0, 


,c:/,+0.891 
■-'-'^^-0.284 


0, 


■'J' "-0.550 


0, 


m 51 +0-4422 
■"-^"-^-0.0095 


3, 


.33 


CASE 


C3 


-19 c + 10.0 


0, 


o„f, + 1.639 
■^'0-1. 135 


0, 


(■9^+3.200 
■'^^'^-1.040 


0, 


ni /tQ+0.5524 
■"-^^''-0.0098 


3 


.26 


CASE 


C4 




0, 


Q 1 Q~t"21 . 888 
■^-^^-0.255 


0, 


f.9^+0.500 


0, 


00261+°-°°^" 
•""^"^-0.00773 





.57 


CASE 


Dl 


or, f, + 160.0 


0, 


■""^-0.059 


0, 


nr,t: + 0.090 
■'^'-'°-0.050 


0, 


onsi+01054 

.UU01_Q QQ27 





.94 


CASE 


D2 


75 n+™-° 

/O.U_gQ g 


0, 


■^"'^-().(i23 


0, 


99C+0.850 
■^^•-'-1.000 


0, 


(,-,9. +0.1679 
■"-^^'-0.0078 


3 


.72 


CASE 


D3 




2, 


g, -,+7.894 
"^^-1.801 


1, 


9oc;+3.780 
■^■^"-1.890 


0, 


(,-, ot;+0.0348 
■"^■^"-0.0090 


4, 


.06 


CASE 


D4 


7 Q+6.0 
'•"-6.0 


0, 


OS4+1.243 
•"^°^-0.512 


0, 


•-^03-0.840 


0, 


(j(,7(,+0.7796 


2 


.10 


CASE 


El 


c + 105.0 
^'■^^-25.0 


0, 


„9r-+7.087 


0, 




0, 


nn35+o-''i93 


1 


.47 


CASE 


E2 


09 c:+45.0 
•^^■'^-25.0 


0, 


49Q+28.250 
■ ^^-0.381 


0, 


Qf;K + 0.540 
■"""-0.420 


0, 


0Q99+0.6392 
■""^^-0.0058 


3 


.85 


CASE 


E3 


10 c: + 12.() 
iO.O_j2() 


1, 


9-, 9+3.718 
■^-^^-0.881 


0, 


f,7. + 1.350 

■"'"-1.700 


0, 


■UUoO-o.oofio 


3 


.63 


CASE 


El 


(i 0- ^" 

.1,(1 


0, 


.,„--|--().771 
il.:!.17 


0, 


- _^ +().4.-)n 

■^ '■' II. -11)0 


0, 




1 


,76 



* We consider both bound and ejected stars. Only unresolved encounters are neglected in this Table. 

^ In units of km s"^. Peak value of the velocity at infinity of the star cluster. The dispersion around the peaJs value is calculated 
considering those values which contain 50% of the total area descending from the peak. 

A,/ = (.//"^ — ./J"), where ./J" and represent respectively the modulus of the initial and the final angular momentum of the 

cluster star. (^^) is the mean value of the variation of the absolute value of the angular momentum of the cluster star, normalized to 

its initial value J^" . The dispersion around the mean value is calculated considering those values -which contain 34% of the total area in 
the left and right wings, respectively. 

is the peak value of . The dispersion around the peak value is calculated considering those values which contain 50% of the 
total area descending from the peak. 

« (^^iSH ) is the mean value of the variation of the binding energy of the binary IMBH, normalized to its initial binding energy SJ,'!,. 

^BH 

The dispersion around the mean -value is calculated considering those values which contain 34% of the total area in the left and right 

wings, respectively. 

* is represents the hardening factor and is given by is = ^ aEbh ^ (cMP 2003). 



the bulk of the supra-thermal stars arc produced. It has 
a non monotonic trend -with the hardness of the binary: a 
harder binary has little excess of angular momentum in its 
initial state relative to that of the incoming star; so the ef- 
ficiency of angular momentum transfer reduces. A loss hard 
binary (having a larger Jbh) is also less efficient, since it 
does not alter much the post-encounter asymptotic velocity. 

We can further quantify the importance of angular mo- 
mentum transfer from the binaxy to the star, by computing 
the fraction of stars which increase J* after an interaction; 
this is given in the last column of Table 4. The fraction of 
stars that do so, over the total (in the sample of the bound 
stars), is significantly high: it is above 70% in many cases. 

Is alignment induced in the scattering process ? In Ta- 
ble 4 we compare the percentage of bound stars which were 
corotating (i.e. for which the scalar product between their 
angular momentum and the angular momentum of the bi- 



nary is positive, that is J, • Jbh > 0) before the 3-body in- 
teraction with the percentage of those stars which are coro- 
tating after the 3-body interaction. Whereas the fraction 
of corotating stars before the interaction is ~ 50% (as one 
can expect given the initial sampling of the data), the frac- 
tion of corotating stars after the interaction is often over 
60%, indicating a tendency of stars to align their angular 
momentum with that of the binary IMBH. This tendency 
is greater if the binary has a big reduced mass (cases with 
A/i = 100, M2 = 50 M0 and with A/i = 50, I'h = 50 Mq) 
and if the binary is moderately hard (i.e. its orbital separa- 
tion is neither too small, because in this case Jbh is compa- 
rable with Jl", nor too largo, since in this case the binary 
would be too soft, and so the kinetic energy exchange would 
be negligible). In case D3 the bound stars which after the 
interaction are corotating with the binary are ~ 70%. This 
means that we can observe a family of supra-thermal stars 
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Table 3. Statistics of the Outgoing States. 



Bound stars (%) Ejections (%) Unresolved Encounters Supra-thermal stars (%)'^ 



CASE 


Al 


45.02 


53.12 


1.86 


21.46 


CASE 


A2 


51.32 


47.10 


1.58 


35.56 


CASE 


A3 


97.60 


1.65 


0.75 


17.23 


CASE 


A4 


99.98 


0.00 


0.02 


0.06 


CASE 


Bl 


53.78 


43.44 


2.78 


13.12 


CASE 


B2 


21.12 


77.90 


0.98 


15.58 


CASE 


B3 


83.36 


16.30 


0.34 


43.36 


CASE 


B4 


99.96 


0.00 


0.04 


0.36 


CASE 


CI 


60.46 


34.82 


4.72 


15.10 


CASE 


C2 


51.30 


45.58 


3.12 


35.24 


CASE 


C3 


95.56 


1.80 


2.64 


17.22 


CASE 


C4 


99.86 


0.00 


0.14 


0.01 


CASE 


Dl 


28.04 


67.84 


4.12 


14.38 


CASE 


D2 


14.44 


84.18 


1.38 


10.68 


CASE 


D3 


75.20 


24.32 


0.48 


44.04 


CASE 


D4 


99.84 


0.04 


0.12 


3.46 


CASE 


El 


27.60 


62.54 


9.86 


17.86 


CASE 


E2 


45.84 


47.42 


6.74 


31.54 


CASE 


E3 


88.68 


2.36 


8.96 


18.42 


CASE 


E4 


98.48 


0.02 


1.50 


0.38 



* We define unresolved interactions those runs where the separation between the star and the binary IMBH never exceeds 30 a after 
10® time steps. 

We define supra-thermal stars those with velocities at infinity comprised between 20 and 40 km s~ ^ . 



Table 4. Statistics for bound stars*. 





Coroti„ (%)'> 


Corot/i„ (%)<= 


More aligned stars (%)'^ 


j/*" > {%)" 


CASE Al 


40 


42 


52 


64 


CASE A2 


47 


59 


61 


72 


CASE A3 


49 


68 


65 


87 


CASE A4 


50 


51 


50 


55 


CASE Bl 


46 


46 


50 


59 


CASE B2 


57 


63 


53 


64 


CASE B3 


50 


58 


56 


83 


CASE B4 


50 


58 


55 


58 


CASE CI 


46 


46 


51 


59 


CASE C2 


45 


49 


58 


64 


CASE C3 


49 


62 


64 


88 


CASE C4 


50 


52 


51 


59 


CASE Dl 


40 


40 


51 


64 


CASE D2 


55 


60 


52 


71 


CASE D3 


51 


70 


62 


95 


CASE D4 


50 


69 


63 


71 


CASE El 


40 


40 


50 


61 


CASE E2 


42 


45 


60 


65 


CASE E3 


48 


53 


61 


92 


CASE E4 


50 


54 


59 


63 



* In this Table we consider only the stars which after the interaction with the binary remain bound to the cluster. 

Percentage (respect to the total of bound stars) of stars which, before the interaction, are corotating with the binary (i.e. for which 
the scalar product between their angular momentum and the angular momentum of the binary is positive). 

Percentage (respect to the total of bound stars) of stars which, after the interaction, are corotating with the binary, independently 
from the initial orientation of their angular momentum. 

Percentage of bound stars which, after the encounter, reduce the angle between their angular momentum and that of the binary. 

Percentage of bound stars which, after the encounter, increase the absolute value of their angular momentum. 



The fingerprint of binary intermediate mass black holes 7 




0.005 



Figure 3. Angular momentum distribution of the cluster star 
in the case D2 {lejt panel), D3 [central panel), D4 (right panel). 
The angular momenta are in dimensionless units (i.e. normal- 
ized to the length, mass and time scales). Dashed line indicates 
the initial angular momentum distribution, solid line indicates 
the post-encounter angular momentum distribution. Evidence of 
angular momentum transfer from the binary to the cluster star 
occurs when the binary has an semi-major axis of 100 AU (case 
D3); while the softest binary (case D4, a=1000 AU) is the least 
efficient from this point of view. On the y-axis the number of 
cases for each bin is normalized to the total number of resolved 
runs (Nfot). 




for 



Figure 4. The histograms show the distribution of -j — 
the case D3 (Mi = 100 M©, M2 = SOM©, a = 100 AU). Solid 
line indicates the distribution after the encounter; dashed line in- 
dicates the distribution before the encounter. The left panel rep- 
resents the distribution of ''^bh j^j. j^ound stars which before 
the encounter were corotating, the central panel the distribution 
for bound stars which before the encounter where counterrotating 
and the right panel the sum of the two distribution, i.e. the dis- 
tribution for all the stars which after the encounter remain bound 
to the cluster. On the y-axis the number of cases for each bin is 
normalized to the total number of bound stars {Abound)- 



which have an anisotropy in their orbital angular momentum 
relative to the center of the cluster, 70% of them rotating in 
one direction and the remaining 30% rotating in the other. 
Obviously, we have to take into account that orbits initially 
are isotropically distributed and that we observe them in 
projection; then the detection of this phenomenon is not so 
immediate. 

Table 4 indirectly shows also another interesting effect. 
The first column indicates that generally less than 50% of 
the initially corotating stars remain bound to the cluster, 
while the post-encounter percentage of corotating stars is 
more than 50%. This means not only that a fraction of ini- 
tially counter-rotating stars becomes corotating, but also 
that an initially corotating star is more easy ejected from the 
cluster than a counter-rotating star. This fact can be intu- 
itively explained considering that, when a counter-rotating 
star interacts with the binary, its relative velocity with re- 
spect to the lighter BH is higher than in the case of a coro- 
tating star, and, then, the cross section is lower. 

To better constrain orbit alignment in the direction of 
Jbh, we considered separately orbits which before the en- 
counter (as initial condition) were corotating (Table 5) and 
orbits which before the encounter were counterrotating (Ta- 
ble 6). For both we investigated the final distributions, and 
in particular we compared the fraction of stars initially coro- 
tating which after the interaction remain corotating (Table 
5, 2nd column) with the fraction of stars initially counterro- 



tating which after the interaction become corotating (Table 
6, 2nd column). We found that most (at least 70%) of the 
initially corotating stars remains corotating, even if during 
very energetic interactions -with the hardest binaries- the 
star tends to forget its initial angular momentum orienta- 
tion (case D2). In contrast, a high percentage of initially 
counterrotating stars becomes corotating, fiipping their an- 
gular momentum. The combination of these two tendencies 
(i.e. the tendency to remain corotating for initially corotat- 
ing stars and the tendency to become corotating for initially 
counterrotating stars) provides the strongest evidence of an 
alignment between and J^" after the interaction (Fig. 
4). In Fig. 5 we give a complete overview of the cases consid- 
ered. 



4 DETECTING SUPRA-THERMAL STARS 

We want now to estimate how many supra-thermal stars are 
produced by interactions with a binary IMBH and if they 
are observable with present instrumentation. 

4.1 Estimating the number of supra-thermal stars 

Since supra-thermal stars have excess kinetic energy and an- 
gular momentum, they are not in thermodynamic equilib- 
rium with respect to the rest of the cluster. Thus, two-body 
relaxation will try to obliterate such anisotropics. Friction 



8 M. Mapelli, M. Colpi, A. Possenti & S. Sigurdsson 

Table 5. Statistics for initially corotating bound stars*. 





Corot/i„ 


Counterrotj^i^ (%)'^ 


More aligned stars (%)'^ 


j/™ > 4" (%)<= 


CASE Al 


93 


7 


49 


70 


CASE A2 


77 


23 


44 


77 


CASE A3 


79 


21 


48 


90 


CASE A4 


82 


18 


42 


53 


CASE Bl 


97 


3 


50 


63 


CASE B2 


67 


33 


34 


67 


CASE B3 


60 


40 


35 


84 


CASE B4 


62 


38 


33 


61 


CASE CI 


98 


2 


50 


62 


CASE C2 


81 


19 


46 


73 


CASE C3 


86 


14 


48 


91 


CASE C4 


86 


14 


45 


60 


CASE Dl 


92 


8 


50 


68 


CASE D2 


65 


35 


31 


67 


CASE D3 


72 


28 


41 


93 


CASE D4 


75 


25 


43 


71 


CASE El 


94 


6 


48 


63 


CASE E2 


86 


14 


48 


73 


CASE E3 


96 


4 


35 


95 


CASE E4 


94 


6 


53 


70 



* In this Table we consider only the stars which before the interaction were corotating. 

Percentage (respect to the total of initially corotating bound stars) of stars which, after the interaction, remain corotating with the 
binary. 

Percentage (respect to the total of initially corotating bound stars) of stars which, after the interaction, become countcrrotating 
respect to the binary. 

Percentage (respect to the total of initially corotating bound stars) of stars which, after the encounter, reduce the angle between their 
angular momentum and that of the binary. 

" Percentage of bound stars (respect to the total of initially corotating bound stars) which, after the encounter, increase the absolute 
value of their angular momentum. 



will return them to the core where they become unrecogniz- 
able, unless some memory of the angular momentum align- 
ment (e.g. Section 3.2) is retained during the relaxation pro- 
cess. 

The existence of supra-thermal stars requires: 

(i) the survival of the IMBH as a binary in the cluster core; 

(ii) the persistence of the suprarthermal stars in the GC 
halo. 

We can investigate these questions through a comparison 
between relevant timcscales. The first timescale to consider 
is the hardening time of the binary IMBH (Quinlan 1996), 
given the stellar number density in the core n and the binary 
separation o 



(4) 



2na^E Gmn 

The second relevant timescale is the relaxation time at the 
half mass radius 



irh — 



ln(0.4iV) \GM 



1/2 



(5) 



where N and M are respectively the number of stars and 
the total mass of the cluster. In Figure 6 we compare the 
two timescales for a = 8.5kms~^, ~ 1 (see Table 2) 
and n = 10® stars pc~*. This plot shows that, when the 
binary is not very hard, i.e., when a > 4 AU, the harden- 



ing time is shorter than the half-mass relaxation time. This 
implies that we can simultaneously observe all the supra- 
thermal stars in the halo for a time comparable to trh with 
no dilution. On the other hand, when the binary has or- 
bital separation a < 4 AU, trh is shorter than thard, and, 
then, we can not simultaneously observe all the produced 
supra-thermal stars in the halo, a part of them being al- 
ready thermalized. However, this effect is balanced by the 
fact that the binary remains in this very hard configuration 
for most of its life, and, therefore, we can observe a consider- 
able fraction (~ 60%) of the total number of supra-thermal 
stars even when the binary is very hard (a < 4 AU) . 
For completeness, we also considered a third timescale: the 
gravitational wave timescale, tg^, i.e. the characteristic time 
for a binary IMBH to coalesce due to gravitational wave 
emission. When this timescale, defined as (Peters 1964; 
Quinlan 1996) 



c^a*(l-e-^) 



2x7/2 



256 G3MiM2(Mi + M2 



(6) 



becomes shorter than thard, the binary would shrink, due 
to gravitational wave emission, at a rate shorter than what 
we expected from thard- Figure 6 shows the timescale tg-u, 
computed for a constant eccentricity e = 0.7 (the initial ec- 
centricity we chose for our runs) and for the range of masses 
considered in this paper. We expect that the assumption of 
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Table 6. Statistics for initially counterrotating bound stars''. 





Corot/i„ 


Counterrotj^i^ (%)'^ 


More aligned stars (%)'^ 


j/™ > Ji" (%)<= 


CASE Al 


7 


93 


54 


60 


CASE A2 


43 


57 


75 


69 


CASE A3 


57 


43 


81 


85 


CASE A4 


21 


79 


58 


57 


CASE Bl 


3 


97 


49 


57 


CASE B2 


57 


43 


79 


60 


CASE B3 


56 


44 


77 


82 


CASE B4 


55 


45 


79 


55 


CASE CI 


34 


66 


52 


57 


CASE C2 


23 


76 


68 


57 


CASE C3 


39 


61 


79 


84 


CASE C4 


17 


83 


57 


58 


CASE Dl 


6 


94 


52 


61 


CASE D2 


54 


46 


78 


76 


CASE D3 


67 


33 


84 


97 


CASE D4 


62 


38 


83 


71 


CASE El 


5 


95 


52 


59 


CASE E2 


15 


85 


69 


59 


CASE E3 


13 


87 


85 


90 


CASE E4 


15 


85 


65 


57 



* In this Table we consider only the stars which before the interaction were counterrotating. 

Percentage (respect to the total of initially counterrotating bound stars) of stars which, after the interaction, become corotating with 
the binary. 

Percentage (respect to the total of initially counterrotating bound stars) of stars which, after the interaction, remain counterrotating 
respect to the binary. 

Percentage (respect to the total of initially counterrotating bound stars) of stars which, after the encounter, reduce the angle between 
their angular momentum and that of the binary. 

Percentage of bound stars (respect to the total of initially counterrotating bound stars) which, after the encounter, increase the 
absolute value of their angular momentum. 



constant eccentricity is correct for star-binary IMBH inter- 
actions, since, when Mi -|- 3> m, the binary eccentricity 
suffers very small changes. We also checked whether this as- 
sumption is correct on the basis of our simulations. We found 
that the mean eccentricity e of the binary IMBH after a 3- 
body encounter is always 0.70, with a very narrow spread 
(the standard deviation a being always 1 — 5 x 10~^), and 
the maximum post-encounter eccentricity is always ^ 0.8, 
in agreement with our statement. 

From Figure 6 we note that tgu, is longer than thard, for 
a ^ 0.5 — 0.3 AU. This implies that the lifetime of a binary 
IMBH is controlled by 3-body scattering on the hardening 
time thard, as long as the binary separation a ~ 0.1 AU. 
Therefore, we can neglect the effect of gravitational wave 
emission, since we consider binaries wider than 0.1 AU. 

The expected number of supra-thermal stars is a frac- 
tion / of the total number N of stars strongly interacting 
with a binary IMBH. The latter is given by (cfr. CMP): 

m^E Vast/ 

where Oo is the initial semi-major axis of the binary and 
Ost is the minimum semi-major axis for the encounter to 
give a supra-thermal star. In eq. (7) we impose ao = 2000 
AU, corresponding to the orbital separation below which the 
binary becomes reasonably hard to generate suprarthermal 
stars, and Ost = 0.1 AU, the orbital separation below which 



the cross section for 3-body encounters becomes negligible. 
We then calculate the number of stars that remain bound 
to the cluster and the number of supra^thermal stars using 
the statistics derived in our simulations (see the percentages 
given in Table 3, respectively in first column -for the bound 
stars- and in fourth column -for the supra-thermal). The 
results are shown in Table 7. The total number of bound 
stars is always greater than 400 and the number of stars 
stirred up to supra-thermal velocities is always greater than 
100, even for the lightest system (Mi= 50 Mq, M2= 10 

Mq). 

Supra-thermal stars may be recognized from their 
proper motion and/or Doppler line shift. Hence we have 
to take into account projection effects of velocity vectors. 
These effects depend in turn on the orientation of the angu- 
lar momentum of the binary IMBH Jbh with respect to the 
line of sight. We find that about 80% of the supra-thermal 
stars can be recognized as such in the best case, i.e. when 
the line-of sight happens to be nearly parallel to the orbital 
angular momentum of the binary IMBH Jbh (if we are mea- 
suring proper motions) or when the line-of-sight is nearly 
perpendicular to Jbh (if we are measuring Doppler-shifts). 
On the contrary, the percentage reduces to about 50% when 
we are in the most unlucky cases, i.e. when the line of sight 
happens to be nearly perpendicular to the orbital angular 
momentum of the binary IMBH Jbh (if measuring proper 
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Table 7. Number of supra-thermal stars^ 



Mb 


Mb 


Interacting stars'^ 


Bound stars'^ 


Supra-thermal stars 


Supra-thermal stars (last 2 Gyrs)° 


50 


10 


728 


467 


128 


86 


50 


50 


1651 


1009 


219 


146 


100 


10 


2060 


1567 


238 


157 


100 


50 


1722 


727 


262 


158 


200 


10 


1945 


1057 


306 


176 



^ These numbers are calculated for binaries having ao= 2000 AU and ast= 0.1 AU 

b In units of solar masses. Mi and M2 represent respectively the mass of the more and of the less massive black hole in the binary 
IMBH. 

Total number of stars which have an interaction with the binary IMBH. 
Number of stars which after the interaction remain bound to the cluster. 
" Supra-thermal stars formed in the last 2 Gyrs. 
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Figure 5. The histograms show the distribution of 7* Jq" f"'^ 
all the considered cases. Solid line indicates the distribution after 
the encounter; dashed line indicates the distribution before the 
encounter. From the up left panel, going toward right: case A2, 
A3, A4, B2, B3, B4, C2, 03, 04, D2, D3, D4, E2, E3, E4. We do 
not plot cases Al, Bl, 01, Dl and El because they are similar 
to cases A2, B2, 02, D2 and E2. On the y-axis the number of 
cases for each bin is normalized to the total number of bound 
stars (Nbound)- 



motions) , or when the line of sight is parallel to Jbh (if mea- 
suring Doppler shifts) . Assuming these percentages, we find 
that the number of recognizable supra-thermal stars over the 
entire cluster is ~ 60-100 for a "light" binary IMBH (A/i = 
50 Mq, M2= 10 Mq) and ~ 130-210 for a "moderately mas- 
sive" binary IMBH (Mi= 100 Mq, Mz^ 50 Mq). 

Our calculation refers to all stars, including also com- 
pact remnants (neutron stars, white dwarfs) and stellar 
types which are too faint for allowing a measurement of 
proper motion or radial velocity. Using the numerical code 
which will be described in the next Section, we find that the 
stars with mass from 0.6 to 0.9 Mq (a mass range which in- 



£,E = Ij ^Tid = 8.5 km s 
Tavitational wave timescale t„ 



Figure 6. Half mass relaxation time, t^j,, and hardening time, 
^hardt function of the initial orbital separation a (in AU) of 
the binary IMBH, for n = lO^pc" 
The shadowed area indicates the ^ 
for binary IMBHs in the considered range of masses (from 60 to 
210 Mq) and with eccentricity e = 0.7. 



eludes red giant (RGB), horizontal branch (HB) and bright 
enough main sequence (MS) stars) represent about the 45% 
of the stars enclosed within^ 0.1 of a GC like NGC 6752. 
This means that, in the case of a binary of (lOOM© , 5OM0), 
we are able to recognize only 50-100 supra-thermal stars. 

There is a further problem concerning the number of 
surviving supra-thermal stars. Nearly all the proposed mech- 
anisms of formation of binary IMBHs in GCs (Miller & 
Hamilton 2002; Sigurdsson & Hernquist 1993) predict that 
such binaries are born within the first Gyr since the forma- 



We consider the region within 0.1 Tc, because we are interested 
in those stars which have the highest probability of interacting 
with the central binary IMBH. 
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tion of the GC itself. If this is true, considering the hardening 
time thard plotted in Figure 6, the binary IMBHs which still 
survive today must have orbital separation ^ 1 AU (for a 
cluster with n = 10^ pc~^). This means that all the supra- 
thermal stars produced when the binary had orbital sepa- 
ration a ^ 3 AU have already been thermalized, because a 
time much longer than trh has elapsed (Fig. 6). Then, we 
can observe only the supra-thermal staxs produced in the 
last ~ 2 Gyr (about 3 trh, because a supra-thermal star can 
have apo-ccntcr distance larger than the half mass radius 
and, therefore, relaxation time longer than trh) - Luckily, the 
orbital separation of a IMBH binary remains in the range 
0.1-3 AU for a long time, because the hardening rate is lower 
for such a hard binary, and a large fraction of supra-thermal 
stars (~ 60%) is produced in this stage. The number of these 
"surviving" supra-thermal stars is reported in Table 7, 6th 
column. Let us consider again the IMBH binary of (100,50) 
A^Q. The total number of supra-thermal stars produced by 
this binary in the last 2 Gyrs is ~ 160. This means that, 
taking into account projection effects and considering only 
bright enough stars (i.e. in the mass range 0.6-0.9 M©), we 
are finally left with only 30-60 recognizable supra-thermal 
stars. 

In this discussion we have not considered the effective 
instrumental errors so far. We now briefly report on them, 
without entering in details. Even if STIS is no more operar 
tive, its accuracy remains a good lower limit for future spec- 
trographs. Observing with STIS stars in a globular cluster 
with distance from the Sun of the order of 4-10 kpc , one can 
expect an error of 1-2 km s"'^ in the determination of radial 
velocity. For example, van der Marel et al (2002), reported 
spectra of about 130 stars in the core of M15 (distance from 
the Sun about 10 kpc) with an observational error of the or- 
der of 1.3 km s~^ . A somewhat higher error can be estimated 
for proper motion measurements with the HST/WFPC2. 
Drukier et al. (2003) combined two sets of observations with 
the WFPC2 (one taken in 1994, the second in 1999) for a 
sample of 1281 stars in NGC 6752, estimating a median er- 
ror of 0.31 mas yr~^ with a mode of 0.17 mas yr~^. For the 
distance of NGC 6752 this means a median error of ~ 6 km 
s~^ with a mode of ~ 3 km s~^, which is still an acceptable 
accuracy to distinguish supra-thermal stars. Similar consid- 
erations hold for HST/ACS (see e.g. Anderson 2002; Ander- 
son & King 2003). Since we defined supra-thermal the stars 
with projected velocity higher than ~ 12 km (16 km 
s~^ for proper motion measurements), an error of 1-2 km 
s^^ (3-6 km s^^ for proper motion measurements) is suffi- 
cient to distinguish supra-thermal from other cluster stars, 
in clusters like NGC 6752. In summary, the main problem 
in detecting supra-thermal stars is not the error on the sin- 
gle measurement but the possibility of observing a sufficient 
large sample of stars, since supra-thermal stars are expected 
to be a very small fraction of cluster stars. 



4.2 Spatial distribution of the supra-thermal stars 

Since supra-thermal stars are only few tens, the possibility 
of recognizing them may be significantly enhanced if their 
radial distribution shows some characteristic feature. Thus, 
it is of interest to study how supra-thermal stars evolve in 
the cluster and what is their radial distribution. 

To this purpose we have explored the dynamics of 




Figure 7. Final projected distribution of the high velocity supra- 
thermal stars {solid line) in the case D3. The peak is at 3 rc- On 

the y-axis the number of cases for each bin is normalized to the 
total number of supra-thermal stars. The (dotted line) shows the 
radial distribution of red giant stars in NGC 6752 (normalized to 
the total number of red giant stars in the sample). 



supra-thermal stars in a cluster under the action of dynam- 
ical friction and the influence of two-body relaxation effects 
using an updated version of the code described in Sigurds- 
son & Phinney 1995. The code first generates a cluster back- 
ground model, i.e. a multi-mass King density profile which in 
our case reproduces that of NGC 6752, and we inject in this 
background supra-thermal stars whose initial positions and 
velocities are those obtained by our 3-body simulations. We 
followed the dynamical evolution of these stars for a random 
time t uniformly distributed in the range < t < 3 trh- In 
this way we can see how supra-thermal distribute in the clus- 
ter before thcrmalization. After evolution for a time t ^ 3 trh 
we select the stars that still are supra-thermal. We applied 
this procedure to the cases Dl, D2, D3 and D4. 

The stars which were still supra-thermal when the sim- 
ulation stopped and that can be observed as supra-thermal, 
taking into account two-dimensional (one-dimensional) pro- 
jection effects and ejections, are about 70% (60%) of the 
initial sample. These stars present a radial distribution like 
that shown in Figure 7 (for the case Dl). This distribution 
is peaked around a radius Vpeak which is of the order of 
few core radii for all the considered systems. For a binary 
IMBH with Ml = 100 M©, M2 = 50 M© and a = 1 AU 
(Dl) rpeak = 3rc and we have calculated that it is enough 
to measure radial velocities within a distance of 6 rc from 
the cluster center, in order to avoid poissonian fluctuations 
to smear out the peak. 

Taking again into account both projection effects and 
luminosity criteria, only ~ 30 supra-thermal stars are lo- 
cated within 6 fc- We found that, if the velocity distribution 
is Maxwellian, ^ 10 stars with velocity in the supra-thermal 
range are predicted to inhabit within a distance of 6 rc from 
the center of NGC 6752. Then, the number of supra^thermal 
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stars produced by interactions with a binary IMBH prevails 
on the high velocity tail of the Maxwellian distribution. On 
the other hand the presence of this tail further dilutes the 
signature of supra-thermal stars. 

In Fig. 7 we have finally compared the radial distri- 
bution of 1984 red giant stars (dotted line) observed in 
NGC 6752 (591 of them from HST/WFPC2 observations, 
1393 from the ESO/MPG WFI. See Sabbi et al. 2004) with 
that calculated for supra-thermal stars (solid line). RGB 
stars appear to be more concentrated toward the center of 
the cluster (the maximum being between and 2 rc); but 
the two distributions are quite similar. 

4.3 Detectability of angular momentum alignment 

Wo now discuss the observability of the signature of the 
angular momentum alignment effect. Column 2 of Table 4 
tells us that in the most favorite case about 70% of the stars 
which remain bound are corotating with the binary IMBH, 
whereas 30% are counterrotating. This means that, if there 
is a binary IMBH of Mi= 100 Mq and M-2= 50 Mq, we 
should have about 90-150 corotating and only 40-60 coun- 
terrotating supra-thermal stars. On the other hand, angular 
momentum alignment effects are visible only for sufficiently 
wide binaries (a > 10 AU). Thus, unless a binary IMBH is 
formed recently in a cluster, the alignment effect today is 
completely washed out by dynamical friction. 

However, we notice that a very low probability mech- 
anism of binary IMBH formation in GCs in recent epochs 
exists. In fact, there is some possibility that a "last single 
BH" (Sigurdsson & Hernquist 1993), ejected from the core 
in the early stages of the GO life, remains in the halo for 
several Gyrs. The relaxation timescale out in the halo is 
long, and if the orbits can circularize there (maybe duo to 
time varying galactic tidal field), then the return time to the 
core is long. When this BH comes back to the core, there 
is a high probability that it forms a binary with the cen- 
tral IMBH. Such a binary IMBH would be originated in late 
epochs. However, the probability for this process to occur is 
low, because it requires the BH to receive just the fine-tuned 
post-encounter velocity needed to remain in the outer halo: 
if the velocity is slightly too high, the BH will bo ejected 
from the entire cluster; whereas, if the velocity is low, the 
return time to the core will be too short. 

The process described here is an unmistakable feature of 
a binary IMBH, which can not be produced by a distribution 
of stellar binaries interacting via 3-body encounters with 
cluster stars. Instead, we need dedicated 3-body simulations 
to check whether even stellar binaries can produce effects 
like those described in §3.1 and §4.2. 



(ii) A fraction of stars tend to align their angular momen- 
tum, J*, with that of the binary IMBH, Jbh- This fraction 
depends on the reduced mass /i and on the semi-major axis 
a of the binary IMBH, but it is always of the order of 55- 
70%. This means that the angular momentum distribution 
of stars that have suffered 3-body interaction with a binary 
IMBH presents a slight anisotropy (§ 3.2). 

(iii) The present theoretical models of binary IMBH for- 
mation indicate that these systems are produced in the first 
Gyr of the GC life. This means that today, in dense clus- 
ters like NGC 6752, binary IMBHs have, due to harden- 
ing, orbital separation a ^ 1 AU. As a consequence, supra- 
thermal stars produced when the binary was wider (more 
than 2 Gyrs ago for the case of NGC 6752) have now al- 
ready thermalized. That reduces the current population of 
residual supra-thermal stars about of a factor two with re- 
spect to the total population produced during the entire 
IMBH lifetime (§ 4.1). On the other hand these "surviving" 
supra-thermal stars do not show any signature of angular 
momentum alignment (§ 4.3), because this effect is evident 
only when the stars interact with a relatively wide binary 
(a > 10 AU). 

(iv) How many supra-thermal stars are expected to be 
produced by a binary IMBH and still surviving in a GC like 
NGC 6752? We estimated few hundreds of supra-thermal 
stars. Once subtracted the fraction of compact remnants 
and faint stars and taking into account projection effects 
and thermalization of the oldest supra-thermal stars (see 
§4.1) we calculate that there may be at most few tens of 
supra-thermal stars suitable to be recognized as such in the 
whole cluster. It can be interesting to notice that there is 
already a claim for the detection of high velocity stars in 
47 Tuc (Meylan, Dubath & Mayor 1991). 

(v) Given the aforementioned relatively small number of 
expected supra-thermal stars it is important to study their 
spatial distribution in order to improve the chances of recog- 
nizing their origin. To this purpose (see § 4.2), we followed 
the dynamical evolution of supra-thermal stars in the clus- 
ter potential (assuming a cluster model which reproduces 
the observational characteristics of NGC 6752), taking into 
account dynamical friction and two body relaxation. We 
found that, before thermalization, suprar-thermal stars tend 
to cluster within 6 Tc- 

(vi) Although we have not performed a detailed simular 
tion of the observability of the few tens of expected supra- 
thermal stars with present detectors, simple considerations 
(see §4.1 and §4.2) suggest that the various observational 
biases (mainly related with the limited field-of-view of the 
most sensitive instruments) may lead to the detection of 
only a sub-sample of the population. 



5 SUMMARY 

In our work we explored the effects that the presence of a 
binary IMBH imprints on the velocity and angular momen- 
tum of cluster stars, duo to 3-body encounters with it. The 
main results of our analysis can be summarized as follows. 

(i) A binary IMBH generates a family of supra-thermal 
stars, i.e. of stars which remain bound to the cluster, but 
have a velocity higher than the dispersion velocity. (§ 3.1). 



In the light of these findings, the search of supra- 
thermal stars and of anisotropics in their angular momen- 
tum distribution can hardly provide useful indications on 
the existence of binary IMBHs in GCs with present in- 
struments. However, it may become a feasible task for fu- 
ture instruments having irmch bigger collecting area and 
equipped with detectors having much larger field of view, 
like (Gilmozzi 2004) the Thirty Meter Telescope (TMT) or 
the Overwhelmingly Large Telescope (OWL). 
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APPENDIX A: CHOICE OF IMPACT 
PARAMETERS 

We chose impact parameters in order to select all the inter- 
actions which have not negligible energetic exchange (about 
ABbh/SITh ~ 10"'^), to allow a complete coverage of all the 
interactions that enter in the supra-thermal domain (20-40 
km s~^). We used the formula of the gravitational focusing 
(Sigurdsson & Phinney 1993), for which 



'2GMtP 



1/2 



(Al) 



where hmax is the maximum impact parameter, Mr the total 
mass of the binary, Vao the relative velocity at infinity and 
p the distance of closest approach. 

We firstly imposed p ~ (10 a) for all the systems (where 
a is the semi-major axis of the binary). Then we made some 
test-run to check what were the best maximum impact pa- 
rameters to include all the strongly interacting stars (and 
then the supra-thermal stars). The results of these checks 
induced us to slightly modify bmax as calculated with eq. 
(Al). In fact, when a binary has a very large semi-major axis 
(case with a=100 or 1000 AU), the choice p ~ (10 a) means 
that the incoming star will pass very far from the center of 
mass of the binary, will interact very weakly with the bi- 
nary (unless it closely approaches the secondary component 
of the binary itself) and never obtain ASbh/Sb"^ > 10~^. 
The velocity at infinity peak is always <^ 10 km/s. 

On the other hand, when the binary has a very small 
semi-major cixis (a < lAU), the choice p ^ (10 a) does not 
take into account all the strong interactions. In fact all the 
stars which have p ~ (10 o) pass so closely to the very hard 
binary to become very fast (> 100 km/s). But, if we consider 
only p ~ (10 a) , we miss most of the interacting stars. In fact, 
if we takep ~ (50—100) o, we still have stars that sufficiently 
approach to the center of mass of the binary to receive an 
outgoing velocity ~ 20-40 km/s (exactly the range of supra- 
thermal stars). And, statistically, stars with p > 10a will 
be much more rmmerous than stars with p < 10 a. Then, 
if we take p ^ (10 a), we miss the bulk of interacting stars 
belonging to the supra-thermal interval. 

Then, both when we choose a small b,nax (p £^ 10 a) for 
a very hard binary and a large bmax (p ^ 10 a) for a wide 
binary, we introduce a bias. In the former case we select only 
high velocity encounters (A£'bh/£'bh ^ 10~^), omitting 
the bulk of interacting stars; in the latter we risk to consider 
also unperturbing flybics (ASbh/^bh < 10"'^). 

Fig. Al represents the post-encounter asymptotic ve- 
locity distribution of the star in the case of a binary with 
Ml =50 Mq, M2=10 Mq, a=100 AU. The dotted line shows 
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Figure Al. Post-encounter asymptotic velocity distribution of 
the star, after the interaction with a binary of Mi =50 Mq, 
i\/2 = 10 Mq, a=100 AU. The solid line shows the case with 
bmax = WO AU (the runs reported in the paper as case A3), 
whereas the dotted line shows a case with 6mai=2000AU (cor- 
responding to adopt p = 10a in eq. (Al)). On the y-axis the 
number of cases for each bin is normalized to the total number of 
resolved runs (Ntot). 
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Figure A2. Post-encounter asymptotic velocity distribution of 
the star after the interaction with a binary of Mi =50 Mq, 
M2=10Mq, a=l AU. The soUd Une shows the case with bmax — 60 
AU (the runs reported in the paper as case Al), whereas the dot- 
ted line shows a case with bmax=i AU (corresponding to impose 
p = la in eq. (Al)). On the y-axis the number of cases for each 
bin is normalized to the total number of resolved runs {Ntot)- 
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the case with bmax='2000 AU (corresponding to adopt p = 
10a in eq. (Al)): most of interactions are very weak fly- 
bies (velocity < 8 km/s), and we cannot count them as 
strong interactions. Instead, the soMd hne shows the case 
with bmax=WO AU (p — a), i.e. the runs reported in the 
paper as case A3: this represents a valid statistic sample of 
all the strong interactions (the pealc being at ~ 13 km/s). 
On the other hand, Fig. A2 shows the post-encounter asymp- 
totic velocity distribution of the star in the case of a bi- 
nary with Mi=50 Mq, M2 = WMq, a=l AU. The dot- 
ted line, showing the case with bmax=l AU (correspond- 
ing p = la), is peaked at very high velocities. Instead, the 
solid line, for bmax=60 AU (case Al), shows that, if we take 
into account also encounters with p > la, but still with 
AEbh/E^ > 10^'^ we have a very different velocity distri- 
bution, peaked at 20-30 km/s. We think that the last one is 
the only statistically significant case, because it takes into 
account all the interacting stars with significant energetic 
exchange. The dotted line, instead, misses the most numer- 
ous class of interactions. 

For these reasons we adopted new criteria for the choice 
of p (to be substituted in eq. Al to derive bmax) as: 

- if a > 100 AU, then p ~ a; 

- if 1 AU< a <100AU, then p ~ (10 a); 

- if a=l AU, then p ~ (60 - 100) a. 

This means, that our selected maximum impact parameter 
is quite constant (~ 100 AU), independent of the semi-major 
axis a of the binary (at least if a < 1000 AU). 

We think that the big maximum impact parameter re- 
quired by systems with a=l AU is due to the uncommonly 
large mass ratio between the binary and the incoming star. 
In fact, when at least one of the two components of the 
binary is so massive, it plays am important role on the grav- 
itational focusing, nearly independently of the semi-major 
axis. Then, the incoming star is attracted toward the center 
of mass of the binary, even if it starts with very large im- 
pact parameter. But a binary of a=l AU has a very small 
geometrical cross-section; then only a very small fraction 
of interacting stars will have an effective pericenter p < a. 
Thus, only a few stars will interact very strongly with the 
binary. Most of the stars will have a pericenter p > a; but 
the binary has a so large binding energy that also these stars 
receive non-negligible post-encounter velocity, and they can 
become supra-thermal. 

Even the behavior of shown in Fig. 2, apparently dif- 
ferent from what predicted by the Heggie's law when a =1 
AU, depends on our choice of bmax- In fact, if for a binary 
with semi-major axis a = 1 AU we choose p ~ a (corre- 
sponding to bmax ^ 1 AU, instead of 60-100 AU), we obtain 
^ 4, higher than the value of derived for wider bina- 
ries and in agreement with Heggie's law. On the other hand, 
we need to select p » a, because we want to include all the 
interactions which have not negligible energetic exchange. 
But this choice lowers the average value of from 4 to 
about 1. 



